5Ж 


УДК 004.89 


СВ. Журба', А.А. Тимченко”, В.В. Скоробрещук', А.М. Гаврилей? 
Черкасский кардиологический центр, г. Черкассы, Украина 
Черкасский государственный технологический университет МОН Украины 


Системный подход 
к диагностированию заболеваний 


В статье рассмотрен процесс научного исследования с точки зрения системного подхода как основного 
инструмента системного анализа сложных объектов и систем в отрасли медицины. Рассмотрены задачи и 
методы диагностирования в задачах управления. Проведена методика системных исследований и 
верификация результатов. В результате математического эксперимента построены прогнозные модели 
составных липидов крови. 


Введение. Основы системного подхода 


Развитие медицины существенно зависит от степени и качества использованных 
интеллектуальных средств. Например, ориентированные на практику математические 
модели, а также математические методы и средства имитации; также компью- 
терные базы данных и знаний, совместно с современными телекоммуникационными гло- 
бальными корпоративными и локальными сетями; использование целевых, организован- 
ных современных информационных средств. 

Процесс научного исследования рассматривается с точки зрения системного 
подхода как основного инструмента системного анализа сложных объектов и систем. 

Стремительное развитие информационных технологий, их использование, в 
задачах наблюдения и прогнозирования требуют разработки и усовершенствования 
существующих методов и создания новых методов исследований и диагностики. 

Диагностика — это специально организованное, систематическое наблюдение за 
состоянием объекта, события, процесса, с целью их идентификации, контроля или 
прогноза. В процессе проведения диагностики могут возникать такие распростра- 
ненные задачи: 

— идентификации параметров моделей; 

— кластеризации данных и результатов; 

— прогнозирование процессов развития. 

Объект исследования можно представить с помощью логической последова- 
тельности: 

<Вход> <><Превращение> <><Выход> 


Исследование решаемости задачи системного проектирования ©) — сложная слабо 
структурированная проблема. Ее нельзя непосредственно декомпозировать на 
определенную совокупность локальных проблем, каждая из которых могла бы быть 
представлена и развязана в узких традиционных постановках инженерного проекти- 
рования. Полученные в процессе исследования решаемости задачи © результаты и 
обнаруженные при этом проблемы относительно исходной информации могут быть 
использованы на ранних этапах проектирования. 

Условия решаемости: 

1) за входом; 
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2) за процессом; 

3) по ресурсу. 

Сложная целеориентированная система: 
1) системные свойства; 

2) этапы жизненного цикла; 


3) ресурсы. 
Логическая схема процесса научного исследования: 


<Системный мониторинг> <><системное проектирование> <> 
<><системные модели (ДДС, КДС, ЛДС)> <> 
<><идентификация (параметрическая и структурная: строение, 
функционирование, развитие)>. 


Процесс научного исследования рассматривается в структуре жизненного 
цикла (ЖЦ) системы как объекта новой техники (ОНТ) [1-3]: 


<системное исследование> <><системное проектирование> <> 
<><изготовление> <><эксплуатация>. 


Объект исследования в процессе системного проектирования и развития должен 
приобретать основные свойства, которые вытекают из концепции 4-х «И»: 


<интеллектуализация> <><информатизация> <><интеграция> <><индивидуализация>. 


Алгоритм системного подхода представляется в виде логической последователь- 
ности [4], [5]: 
<ицель> <><задачи> <><методы> <><алгоритмы> <><средства>. 


Естественными объектами исследования выступают: процессы, явления, события 
или системы реального мира, математическими моделями которых являются диф- 
ференцированные динамические системы (ДДС) [5]. 

Поиск конкретных моделей объектов исследования выполняется методами 
идентификации (структурной и параметрической) [3]. 

Основная технологическая последовательность на основе методов системного 
проектирования есть [4], [5]: <синтез> <><анализ> <><принятие решений>. 


Применение системного подхода 
1 Задача и методика системных исследований. 
Верификация результатов 


В основу методики обработки данных положены методы параметрической и 
структурной идентификации. 

1 Этап структурной идентификации. В основе рассмотрения задачи струк- 
турной идентификации на базе методов структурного синтеза положено определение 
системы как сети связанных элементов. Тогда, соответственно, совокупность связей 
определяется как структура системы. 

Известно несколько подходов к решению этого класса задач; общая постановка 
выплывает из формального определения системы в эволюционном развитии. Методы 
и их использование дают возможность рассматривать задачи структурного синтеза 
современных систем управления как теоретическую основу для последующего 
конструирования и реализации на компьютере. 

1.2 Этап параметрической идентификации. В основе построения моделей 
современных объектов автоматизации чаще всего лежит использование понятия 
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«функция» (в самом широком понимании этого слова). Это базируется на определении 
функции как закона, по которому к значению аргумента ставится в соответствие 
значение функции. Применяют такие методы и алгоритмы, которые составляют основу 
информационной технологии построения моделей (поиску законов функционирования 
в самом широком понимании этого слова): метод избранных точек, метод регрес- 
сионного анализа (наименьших квадратов), алгоритм непрерывного определения 
статической характеристики объекта управления (фильтр Габора), метод экстраполяции 
(прогнозирование), метод группового учета аргументов (МГУА), для построения 
сложных моделей и систем и другие. 

1.3 Этап группового учета аргументов (МГУА) [6], [7]. Свое применение МГУА 
нашел в разных областях знаний, включая структурную, параметрическую идентифика- 
цию и прогнозирование. МГУА базируется на принципах самоорганизации моделей. 
При самоорганизации моделей исследователь задает большое число моделей-пре- 
тендентов и с помощью отдельных критериев определяет модель оптимальной слож- 
ности (МОС). Параметры модели-претендента рассчитываются с помощью хорошо 
известного метода наименьших квадратов на базе информации, которая получена в 
результате наблюдений над объектом исследования. 

Идея метода самоорганизации заключается в проверке полученной модели на 
информации об объекте исследования, что не использовалась при определении ее 
параметров, то есть используется принцип внешнего дополнения. 

1.4 Этап кластеризации [7]. Различают такие большие группы методов кла- 
стеризации: методы иерархической кластеризации, процедуры типа упорядочивания 
матрицы расстояний, процедуры эталонного типа, алгоритмы типа разрезания графа, 
комбинированные методы и другие. Кроме того, методы кластеризации можно клас- 
сифицировать в зависимости от вида данных: кластеризации данных, значения которых 
заданы четко, и данных, значения которых заданные не четко, например, в виде функции 
принадлежности. Также методы кластеризации можно разделить в зависимости от вида 
кластерной структуры: кластеры не пересекаются; кластеры пересекаются; объекты 
образуют иерархическое дерево; отношение произвольной структуры. 


2 Задача и методы диагностирования в задачах управления 
процессом лечения 


Рассмотрим такие задачи диагностики: задачу идентификации параметров 
модели; задачу кластеризации данных; задачу прогнозирования процессов развития. 
Задача диагностирования как задача идентификации. Объект управления 
находится под воздействием нескольких входных величин. Среди них различают: 


вектор регулирующих влияний д=(д,,4>,...М„), вектор возмущающих влияний 
й=(1,,4>,...^„) и А, - оператор функционирования объекта. 


На выходе объекта есть вектор ф = (ф,,$......ф,), одна из компонент которого 


является показателем качества, который, в свою очередь, может быть функцией 
нескольких переменных. 

Статистической характеристикой объекта является зависимость значения 
качества (Ф) от всех входных величин в установившемся режиме 


Ф= (и, А). (1) 
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Уравнения динамики объекта записывают в виде 
Е(РФ=Л( и, А), (2) 


где (р) — полином от р = в. 
& 

Нахождение /(р) является задачей определения /(,А) — задачей определения 
статических характеристик объекта и является задачей идентификации. 

Процесс идентификации можно представить с помощью следующей блок- 
схемы (рис. 1): 


лдс 


ый Ц 
——> Система 2, - 
объект исследования 


Желаемая 
система 


обработки: СКО -— срехнеквапратичная омиока 


Критерии кк - ‘коэффициент корреляции 
ранжировка 


Гипотеза по 
одн. 


х 
Мотоды 
7 идентификации 


Рисунок 1 — Блок-схема алгоритма (процесса) идентификации 


Идентификация параметров моделей сложных систем (ШФЕМТ (РУМ$)) осно- 
вывается на методе наименьших квадратов и позволяет определять параметры моделей 
статических и динамических систем исходя из критерия минимума среднеквадратической 
погрешности. 


Статическая система задается системой функций многих переменных, которые 
подаются в виде степенных разложений к 9-й степени включительно: 


ЕО) = Л(хь хо... хр, Р=Ь..,И. (3) 


В качестве примера статическую характеристику системы с одним возмущающим 
влиянием и одним регулирующим влиянием можно представить в виде кривой 
второго порядка: 


Ф=а, чала, и+а, А? ча. и? +а, ди (4) 
Динамическая система задается линейной системой дифференциальных уравнений: 
Х=АхХ+ВОи. (5) 


Процессы прогнозирования реализуются на базе прогнозных моделей, которые 
строятся по методу Пауэлла (РКОСРОУ): 


Хи = 40 +а.Хи- 42›Х+...НаХ и р (6) 
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Работа из ППШ возможна как в диалоговом автономном режиме, так и с ис- 
пользованием среды МАТЕАВ [4]. 

Заметим, что если сумма квадратов отклонений экспериментальных значений 
от рассчитанных будет больше некоторого заранее заданного числа, то порядок уравне- 
ния необходимо повысить. Задача определения коэффициентов уравнения является зада- 
чей параметрической оптимизации. Если необходимо определить вид уравнения, то такая 
идентификация является структурной. Определение коэффициентов уравнения вычисля- 
ется с помощью метода наименьших квадратов, то есть путем минимизации функции 


1 
2 
О= ре =. ф”“" ) | (7) 
11 
где ф’““Р — рассчитанные значения функции Фф , 


ф7“" — экспериментальные (априорные) значения, 1 = 1,[. 


3 Задача кластеризации 


Дана совокупность наблюдений (объектов, процессов, экспериментов), 


Х=(Х,Х,,..,Х„), т -— количество наблюдений, Х, ={х,}, /=1п, п — количество 


характеристик модели. Все наблюдения относятся к определенной общей категории, 

в рамках которой необходимо выделить (если это возможно) подгруппы (кластеры). 
Пусть К - количество кластеров, на которые необходимо разбить совокупность 

наблюдений (К <<т), С=(С,С.,....С к) — совокупность кластеров (множественные 


числа наблюдений, которые не пересекаются), С; = {с,}, /=1р,, р; - количество 


К 
наблюдений в 1-м кластере У р.=т. 
1=1 
Каждый кластер С; имеет центр (ядро кластера) у,, множество центров 
кластеров обозначим У = {у}, 1=1,К . В пространстве наблюдений определим меру 
их близости (критерий) 4(Х,,Х’,). В качестве критерия может выступать евклидово 
расстояние, корреляция и другие. Необходимо решить две основных задачи: 
К 
1) распределить совокупность наблюдений между кластерами, то есть Х = (с, : 
1=1 


2) определить центры кластеров таким образом, чтобы минимизировать функционал [7]: 
КР: 


У У а(с,»у,) > ша. (8) 


11 75 


4 Метод экстраполяции (прогнозирование процессов лечения) 
Содержание метода экстраполяции заключается в следующем. В первую очередь 

строится класс экстраполяционных функций Г(,а1,а2,....аи), который зависит от 

параметров а,,а,.....а„. Выбор значений этих параметров, которые задают нужную 


экстраполяционную функцию, происходит из условия минимальности суммы квадратов 
разниц значений функций / и Х при заданных значениях #=4,6,,....й,. Указанная 


сумма имеет вид 


5=[/@ ‚а, а,,...а„)-Х@ + @, а ,а,...и„)-Х| +... + аа, а,,...й,)-Х@)Р (9) 
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Поиск значений параметров а,,а......а„, которые превращают величину 5 в 


т? 
минимум, происходит любым из способов нахождения экстремумов. Для ряда простых 
функций (а ,а,,...а„) необходимый минимум может быть найден классическим 


путем в результате решения системы уравнений 
Е М (= 1,2....,т). (10) 


Метод наименьших квадратов (МНК) удобен тем, что позволяет применять 
аналитические приемы нахождения минимума. 
Процессы прогнозирования реализуются на базе прогнозных моделей, которые 
строятся по методу Пауэлла [8], [9]: 
Хы = 40 + Хи + а ХА, +...+а, Хи. (11) 


5 Прогноз содержания липидов в сыворотке крови 

В рамках данной постановки задачи и формулировки цели исследования за- 
ключается в построении прогноза содержания липидов в сыворотке крови на основе 
следующих данных. 

Выходные данные. Выходными данными для построения прогноза является общее 
количество липидов в сыворотке крови, содержание холестерина, фосфолипидов, |- 
липопротеидов, обмен которых изучался в динамике экспериментального инфаркта 
миокарда на 1, 2, 3, 5, 10, 20 сутки (табл. 1). 


Таблица 1 — Начальные данные 


Сутки Общее кол-во, | Холестерин, | Фосфолипиды, | В- Липопротеиди, 
мг % мг % мг % мг % 
норма 522 120 220 160,5 
1 480 170 190 190 
2 480 138 170 200 
3 440 92 150,3 220 
5 448 78,6 180 240 
10 465 83 180 260 
15 466 96 182 231 
20 467 108 185 212 
Используя программную систему Ргогро\ и выражение (11), найдем резуль- 


таты прогнозирования на 25, 30, 35, 40 сутки. 
Данные, полученные в результате эксперимента, представлены в табл. 2. 


Таблица 2 — Результаты эксперимента 


Сутки Общее кол-во, | Холестерин, | Фосфолипиды, | В- Липопротеиди, 

мг % мг % мг % мг % 
25 479.22 99.91 182.21 193.30 
30 486.57 105.02 182.78 178.62 
35 493.70 109.34 183.30 166.50 
40 500.62 112.99 183.75 156.50 
45 507.33 116.06 184.15 148.25 
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Построим полученные прогнозы в графических плоскостях (рис. 2): 


190 


5Ж 
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70 
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Рисунок 2 — Графики прогнозирования содержания липидов в сыворотке крови 
Заключение 


Таким образом, изложена системная методика последовательного анализа развития 
заболевания сердечно-сосудистой системы на примере данных липидного обмена. 
Основу математической обработки составляют экспериментальные данные, которые 
характеризуют содержание липидов в крови на разных стадиях инфаркта миокарда. 
В результате математического эксперимента построены прогнозные модели составных 
липидов крови. Следует отметить, что программа Ргогро\" хорошо работает с данными, 
которые описываются линейной зависимостью, а в нашем случае такой зависимости 
не наблюдалось. 

Представленные прогнозы свидетельствуют о возможности использования анали- 
тических платформ, которые реализуют задачи диагностики в задачах медицины. 
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С.В. Журба, А.А. Тимченко, В.В. Скоробрещук, А.М. Гаврилей 

Системний шдхщ до дагностування в задачах медицини 

У статт! розглянутий процес наукового досл1лдження з точки зору системного шдходу як основного 
1нструмента системного анал!зу складних об’ектв 1 систем в галуз! медицини. Розглянут! задач! та 
методи дагностування в задачах управлёння. Проведена методика системних дослджень та верификаци 
результатв. У результат! математичного експерименту побудован! прогнозн! модели складових липдав крови. 


5.7. Яригра, А.А. Тутсйепнко, И.И. 5Когогехйсйик, А.М. Саугйеу 

Зучет АрргоасВ №0 О1аспозто РгоШет$ ш Медсше 

ш бе агафе Фе ргосез$ оЁ зс1епиЙс гезеагсВ ш {егтз$ оЁР зубет арргоасВ аз ше Баз1с 1ю0] оЁ зузет 
апа[уз1$ оР сотр[ех об]ес6 апа зузетз ш те4!са| зслепсе. ТВе ргоМетз ап4 те о45 ш Фазпозте сопго] 
ргоегл$ аге сопз14егеа. ТЬе плео4 оЁ зужетайс гезеагсЬ ап4 уемйсаНоп оЁ гези$ аге сопдисе4. Аз а 
тези оЁР тафетайса] ехрегитепе ргед@сйуе по4е|5 оЁ Пр! сотропеп{5 оРГЫоо4 Вауе Бееп бий. 


Статья поступила в редакцию 16.07.2011. 
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